On chemical bonding of Helium with hep-Beryllium 
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Chemical inertness is the key property of helium determining its solubility, distribution and accu- 
mulation kinetics in metals. Against all expectations, our ab initio calculations show a substantial 
chemical bonding between He and Be atoms in the hep-Be matrix when He occupies a non-symmetric 
position in a basal plane. 
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Helium is the most inert element of the periodic table. 
Having its hrst electronic level (ls-shell) completely occu- 
pied, the electro-neutral state of the atom is very stable. 
In an excited state, when one of the electrons moves to 
the next energy level, He may form short-living dimers 
with fluorine or chlorine: He-F and He-Cl. The very ex- 
istence of these excimers can be viewed as a proof of the 
degree of chemical inertia of helium. Chemical inertia 
determines the helium behavior in the irradiated struc- 
tural materials of nuclear reactors, where He is formed in 
substantial amounts due to nuclear fission reactions. Ac- 
cumulation of He in metal alloys in the form of the inert 
lattice gas, leads to their so-called "helium swelling" and 
embrittlement. 

We studied the distribution of interstitial He in hexag- 
onal close-packed (hep) beryllium by means of ab initio 
modeling. Employing a pseudopotential technique^, we 
obtained a surprising result: interstitial He forms a chem- 
ical bonding with Be matrix 2 . 

Further, using a highly accurate all-electron full- 
potential linearized augmented plane wave (FLAPW) 
method, we confirmed results of pseudopotential calcula- 
tions, analysed partial DOSes and their integral quanti- 
ties. We performed analysis of chemical bonding in the 
Be- He system using Crystal Orbital Hamilton Population 
(COHP) indicator 3 , implemented in SIESTA package^. 
In this Letter the results of our investigations are pre- 
sented. 

The Density Functional Theory (DFT) was used in 
Generalized Gradient Approximation (GGA£), as im- 
plemented in the Quantum-ESPRESSO (QE) program 
package-. The same GGA approximation was used for 
calculations using FLAPW method (Wien2k software 
package^) . 

We have chosen 22 positions of interstitial He in hep- 
Be (space group 194, two Be atoms in the unit cell at 
positions (1/3,2/3,1/4), (2/3,1/3,3/4)) as starting geom- 
etry configurations for the calculations. These positions 
included the centers of octahedra and tetrahedra, as well 
as the midpoints between the centers of these polyhedra 
(Fig. [TJ. The midpoints between the two neighboring 
lattice sites were also included as testing positions for in- 



terstitial He. All positions were generated by relation: 
1/2 {pi,P2,Pz) + 1/2 ( r i + r j): where summation i,j 
is over octahedra and tetrahedra, and Pi,P2, P3 take val- 
ues {0, 1} independently. 

Calculations were performed for Be64He (4x4x2 hep Be 
unit cells), BeggHe (4x4x3 hep Be unit cells) supercells 
and included relaxation of atomic positions. Atomic nu- 
clei were moved to their equilibrium positions while forces 
acting on nuclei were greater than 1 mRy/Bohr. Lat- 
tice parameters of supercells were according to hcp-Bc, 
that were obtained by minimisation of the total energy 
with respect to a and c/a. Calculated lattice parameters 
from pseudopotential method a = 2.246 A, c/a = 1.573 
are comparable with experimental values a = 2.2858 A, 
c/a = 1.568 7 , whereas FLAPW data a = 2.272 A, 
c/a = 1.565 are in adequate agreement. 

On the basis of the total-energy calculations, we found 
five equilibrium positions of He atom (zero force on He 
nucleus) in the beryllium lattice. Fig.[T]shows these equi- 
librium positions as green and orange circles. To denote 
positions we are using Wyckoff notations^ and fractional 
coordinates (Table fl]). The other 17 studied positions of 
He atom in hep-Be matrix, including the tetrahedral one, 
were found to be non-equilibrium. Among 5 equilibrium 
positions, three positions (b),(d),(ho) are stable (Fig.[]Ji, 
green circles), and two (a),(g) are unstable (saddle points 
with negative curvature of potential energy surface for 
He displacement along [0,0,1]; Fig. QJi, orange circles). 
Both FLAPW and pseudopotentail methods produced 
close values of He-Be bond length (Table . The energy 
of isolated helium atom, which is needed to calculate the 
solution energy of He in Be lattice, was obtained by per- 
forming total energy calculations for simple cubic lattice 
of He atoms with lattice parameter 10 A. Calculated so- 
lution energies of interstitial He are presented in Table HIl 
The previous ab initio results of Ganchenkova et al£ and 
Zhang et aZi 10 ' 11 are also shown for comparison. 

According to all studies^— He solution energy is in 
the range of 5.5-6.5 eV. He solution energy for octahedral 
cavity {(a) position) has the value of 6.0-6.5 eV and is in 
good agreement with the value of 6.1 eV fro m 10 ' 11 and 
close to 6.8 eV reported by Cayphasi 2 - from EAM poten- 
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tial calculations. Solution energies for positions (b),(d) 
are in well agreement between researches. According to 
our present findings, the most energy favorable He posi- 
tion is (ho ) with the solution energy of 5.3-5.8 eV. It was 
obtained by moving He from points h(x — 1/2, z = 1/4/ 
and k(x = 1/2, z = 5/8/ (Fig. [1^, grey circles). In (ho) 
He is located in the basal plane near a strongly displaced 
Be atom (Fig. [TJj) . It is worth to be noted, that this po- 
sition (ho) was never considered in previous studies^— 
as a place for He. 

Zero point vibration energy (ZPE) is not small for in- 
terstitial He. To estimate this quantity we calculated 
phonon frequencies at T-point for BeggHe and hcp-Beg6 
supercells within DFPT approach as implemented in QE 
software package. ZPE was calculated as a sum of delta 
functions, that typically produces 5% uncertainty of ex- 
act value for supercells of such size (see e.g. ref. LL3|). For 
positions (b),(d),(ho) ZPE contributions to the solution 
energy are 0.09, 0.08 and 0.15 eV, that does not change 
energy landscape significantly. This correlates well with 
result by Zhang et al^Q, where ZPE correction was found 
to be in narrow range of 0.12-0.13 eV for all studied po- 
sitions. 

A deep potential well for the (ho ) position, the presence 
of only one Be atom in the first coordination shell of He 
and small He-Be distance (1.29A) for this configuration 
demonstrate the existence of chemical bonding between 
He and Be atoms. To get further insight and investi- 
gate the nature of this bonding, we proceeded with the 
highly-accurate FLAPW method to calculate the density 
of electronic states (DOS). 

Two configurations of He in Be matrix were chosen for 
comparative DOS analysis, (a) and (ho) (Fig. [TJ). Fig- 
ures [2j) and [2fc show partial DOSs of the He-Be system 
for the cases, where He occupies the octahedral cavity 
(a) and the position in the basal plane near the lattice 
site (ho). Partial DOSs for the pure hep-Be system are 
also shown in Fig. [2^l for comparison. It is seen that the 
valence band of He-Be system is formed by the He 2s 




FIG. 1. (Color online) Equilibrium positions of He in hep-Be 
and atomic relaxation near He atom in position (ho ). Stable 
positions (b,d,h®) are shown as green circles, unstable (a,g) 
as orange ones. 
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FIG. 2. (Color online) Partial DOSes of the He-Be system: 
a) pure hep-Be; b) He occupies the octahedral cavity (a); c) 
He at position (ho). 



and the Be 2s- and 2p-states. The occupied He ls-state 
is located 18 eV below the Fermi level. The densities of 
electronic states for He in one of the (a) or (ho) positions 
appear to be substantially different, which reflects the 
different degree of overlap of He electronic charge with 
that of neighboring Be atoms. 

Together with substantial decrease of the He-Be dis- 
tance, the (ho ) configuration shows the increase of hy- 
bridisation of the He s-states with the 2s- and 2p-states 
of Be. The hybridisation of low-lying He ls-state and 2p- 
states of Be also takes place. We need to mention here 



TABLE I. Equilibrium positions of He atom in hep-Be (space 
group 194), number of nearest Be atoms z, distance He-Be 
d(A), calculated by plane-wave pseudopotential (PS, QE) and 
FLAPW (Wien2k) methods for Be^He, BegeHe structures. 
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subscript "0" indicates minimal solution energy among all 
positions 
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TABLE II. Solution energies of helium (eV) in hep-Be from 
DFT(GGA) calculations. VASP - Vienna Ab initio Simu- 
lation Packag e 16 ' 17 ; "Y/N" symbols indicate the presence or 
absence of ZPE contribution. 
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a With PBE 5 exchange-correlation. 

b With PW91— — exchange-correlation. 



that the valence He s-states originate from the empty 
2s° atomic level of helium. When the Be-He system is 
formed, this empty state moves below the Fermi level, 
and becomes partially occupied, forming the s-band of 
Be-He system. These results reveal an increasing hy- 
bridisation of the He electronic states with those of neigh- 
bouring Be atom when helium is placed in (ho ) position. 
Due to the formation of the chemical Be-He bond, when 
the He atom is located in the stable equilibrium (ho) 
position, attraction occurs between He and the near- 
est Be atom. As a result, the He and Be atoms are 
displaced from the centers of local symmetry and the 
distance between them is substantially reduced. The 
strength of the chemical bond is nearly equal to the total 
energy difference between (a) and (ho) configurations, 

Echem ~ e o — e h a ~ 0-7eV. 

The analysis of the partial density of states is not 
enough for complete understanding of the nature of the 
chemical bonding between He atom and Be matrix. We 
performed additional analysis of the chemical bond for- 
mation using the concept of bonding and antibonding 
states within the framework of COHP chemical bond 
indicator—. We used pseudopotential method and local 
atomic orbital basis set, as realized in SIESTA program 
package^. The calculations were done in GGA approxi- 
mation for exchange-correlation functional 5 , and double- 
zeta polarized basis set was used, which included doubled 
2s and 2p orbitals plus polarizing 3d orbital for Be atom, 
and doubled Is plus polarizing 2p orbital for He. An op- 
timization of atomic positions of BeggHe supercell with 
He in positions (a) or (ho) was performed. The obtained 
relaxed positions of Be and He atoms were very close 
to those, obtained by the FLAPW and pseudopotential 
(QE) methods. The calculated Be-He bond lengths for 
positions (a),(b),(ho) were 1.809 A, 1.580 A, 1.276 A, 
which agrees well with the results, obtained by other 
methods (Table [[}. 

We plot in Figure [3] the calculated COHP curves for 
He and Be atoms for He at positions (a), (b) and (ho). 
Figure [3^i) shows COHP curves for He atom and near- 
est beryllium Befjc- The plot demonstrates that only 
bonding states are present for He atom located in (ho ) 
position, which is not the case when He is located in (a) 
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FIG. 3. (Color online) COHP curves for Be96He structure 
with helium at positions (a),(b) and (ho): a) He and nearest 
BeH C ; b) He and beryllium matrix; c) BeHe and beryllium 
matrix 



and (b). Antibonding states near Fermi level originate 
from atomic He 2s-states (see Fig. [2Jd). Integrated val- 
ues (to Fermi level) of these COHP curves for (a),(b),(ho) 
positions are 0.32, -0.19, -3.42 eV. 

COHP curves for He and whole beryllium matrix are 
shown in Figure [3|d. The bonding strength of He ls- 
states is the same for all 3 positions, but reduction of 
He-Bene distance (especially strong for (ho)) decreases 
antibonding states near Fermi level. Integrated COHP 
values are 1.93, 1.38, -1.09 eV, confirming He-Be bonding 
in position (ho ). This fact, as well as the small He-Bene 
interatomic distance, allows us to conclude that a distinct 
molecular Hc-Bcho pah' is formed in Be-hcp matrix. 

Soluted in hep beryllium, helium also changes Be-Be 
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FIG. 4. (Color online) Distributions of electron density 
(e/(a.u.) 3 ) for structure Bee4He with He atom: a) in posi- 
tion (a) and b) in position (ho). 



bonding. Figure [3}: shows COHP curves for Bene atom 
and beryllium matrix, as well as COHP curve for pure 
hep-Be. It can be seen that the curves for (a), (b) posi- 
tions are close to that of hep-Be. For (ho) position, He 
atom induces antibonding Be-Be states near Fermi level, 
thus weakening Be-Be bonding. 

To visualize the chemical bond in Be-He system, we an- 
alyzed the spatial distribution of valence electronic den- 
sity in model Be64He supercell, when He is located in 
(a) or (ho) positions. Figure 2] shows these distributions 
for density values larger than 0.01 e/(a.u.) 3 . When He is 
located in the octahedral cavity (FigHk), we do not ob- 
serve electronic charge concentration between He and Be 
atoms (marked as red and blue spheres respectively) , and 
the electronic density (marked as a green cloud) appears 
to be "pushed" out of the Be octahedron. Completely 
different charge distribution is observed in case of He 
located in (ho) position (Fig. In this case the elec- 
tronic charge "connects" He and Be atoms, which shows 
a covalent character of chemical bonding. 

In conclusion, our ab initio calculations show a sub- 
stantial chemical bonding between He and Be atoms, 
when He occupies a non-symmetric position in a basal 
plane in the hep-Be matrix. This fact should be taken 
into account, while analysing the kinetics of the process 
of helium accumulation in hep-Be. The experimental evi- 
dence of the existence of the Be-He bond may be obtained 
by means of different methods, including the diffraction 
experiments. We believe that our results should stimu- 
late the experimental efforts in this direction. 
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